Collective oscillations of a Fermi gas in the unitarity limit: 
Temperature effects and the role of pair correlations 
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We present detailed measurements of the frequency and damping of three different collective 
modes in an ultracold trapped Fermi gas of ®Li atoms with resonantly tuned interactions. The 
measurements are carried out over a wide range of temperatures. We focus on the unitarity limit, 
where the scattering length is much greater than all other relevant length scales. The results are 
compared to theoretical calculations that take into account Pauli blocking and pair correlations in 
the normal state above the critical temperature for superfluidity. We show that these two effects 
nearly compensate each other and the behavior of the gas is close to the one of a classical gas. 

PACS numbers: 34.50.-s, 05.30.Fk, 39.25.+k, 32.80.Pj 
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I. INTRODUCTION 

The study of collective oscillations in quantum liquids 
and gases has yielded a wealth of insights into the prop- 
erties of strongly correlated systems. An early example 
concerning strongly correlated Fermions is the observed 
transition from ordinary first sound to zero sound in the 
normal state of liquid ^He as the temperature is low- 
ered [H. In this Article we explore related phenomena 
in an ultracold quantum gas of fermions in the unitarity 
limit [2] by measuring three different collective modes 
under similar conditions. The frequency and damping 
of the modes exhibit the characteristic transition from 
hydrodynamic behavior at low temperature to collision- 
less behavior at higher temperature. The experimental 
observations are compared to theoretical model calcula- 
tions that apply to the normal state of the gas above the 
critical temperature Tc for superfluidity. In the unitar- 
ity limit, the strongly correlated normal state between 
Tc and the Fermi temperature Tp is arguably not as well 
understood as the T = superfluid phase Q . It is shown 
that the combined effects of temperature and pair corre- 
lations account for most of the observed features in this 
interesting temperature regime. 

Our measurements of the collective modes are car- 
ried out for an elongated trap geometry, which has pre- 
viously been shown to be well suited for studying the 
dynamical behavior of a strongly interacting Fermi gas 
0, H, [1, 0j B B US] ■ We focus on two collective excitations 
of a cylindrically symmetric cigar-shaped cloud, namely 
the radial compression mode and the radial quadrupole 
mode. In addition we study the scissors mode under con- 
ditions where the cloud exhibits pronounced ellipticity 
in the plane perpendicular to the direction of the cigar- 
shaped cloud. In all three modes, the cloud oscillates 



mainly in the plane normal to the direction of the cigar- 
shaped cloud. For a sketch of the modes see Figure [T] 

Previous experiments on collective modes in a strongly 
interacting Fermi gas studied the effect of the interaction 
strength in the zero temperature limit. ^^^5, 6, 8, 9]. Sys- 
tematic investigations were performed studying the ra- 
dial compression mode [3, lEI, and the radial quadrupole 
mode 9] . Measurements on the compression mode served 
as a sensitive probe for the equation of state of the gas 
in the zero temperature limit throughout the BEC-BCS 
crossover regime. In contrast to the compression mode, 
the frequency of the radial quadrupole mode allows one 
to test the hydrodynamic behavior without being influ- 
enced by the equation of state. This made it possible 
to investigate the transition from hydrodynamic to col- 
lisionless behavior with decreasing coupling strength of 
the atom pairs on the ECS side of the crossover. 

While the hydrodynamic behavior in the zero- 
temperature limit is now well understood as a result of 
superfluidity, an understanding of the effects of temper- 
ature on the collective modes has remained a challenge. 
Only few experiments have so far addressed this problem 
0, 0, B [l3| ■ Previously, the temperature dependence of 
the radial compression mode [tI] and the scissors mode 
was studied [lO|. Our present experiments aim at ad- 
dressing the open questions raised by the different results 
obtained in these experiments: The frequency and damp- 
ing of the radial compression mode was studied as func- 
tion of the temperature in an experiment performed at 
Duke University [7] . There the mode frequency appeared 
to stay close to the hydrodynamic value even for tem- 
peratures exceeding the Fermi temperature. This sur- 
prising finding stands in contrast to scissors mode mea- 
surements, performed later at Innsbruck University p^ . 
which clearly showed a transition to coUisionless behavior 
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Figure 1: Sketch of the three collective modes investigated 
in this work: the compression mode, the quadrupole mode 
and the scissors mode (from left to right). The oscillations 
take place in the plane of tight confinement, perpendicular 
to the direction of the elongated, cigar-shaped cloud. While 
the compression mode represents an oscillation of the overall 
cloud volume, the other two modes only involve surface defor- 
mations. Exciting the quadrupole mode leads to an oscillating 
elliptic shape. The scissors mode appears as an angular oscil- 
lation of an elliptic cloud about a principal axis of an elliptic 
trap geometry. 

in the same temperature range. Furthermore the Duke 
data on the damping of the compression mode did not 
show a maximum as it was seen in the Innsbruck data on 
the scissors mode measurement. These apparent discrep- 
ancies are a particular motivation for our present study 
of different coUective modes under similar experimental 
conditions. 



II. EXPERIMENTAL PROCEDURE 

The apparatus and the basic preparation methods for 
experiments with a strongly interacting Fermi gas of ^Li 
atoms have been described in our previous work pU . [l^ . 
As a starting point, we produce a molecular BEC of ^Li2. 
By changing an external magnetic field, we can control 
the interparticle interactions in the vicinity of a Feshbach 
resonance, which is centered at 834 G The measure- 
ments of the collective modes are performed at the cen- 
ter of the Feshbach resonance, where the interactions are 
unitarity limited. 

The atoms are confined in an elongated, nearly har- 
monic trapping potential, where the trap frequencies lOx 
and LOy in the transverse direction are much larger than 
the axial trap frequency lOz- The confinement in the 
transverse direction is created by an optical dipole trap 
using a focused 1030 nm laser beam with a waist of 47 /im. 
Note that the Gaussian shape of the laser beam leads 
to significant anharmonicities in the trapping potential. 
The potential in the axial direction consists of a com- 
bination of optical and magnetic confinement; the mag- 
netic confinement is dominant under the conditions of 
the present experiments. The trap parameters, given 
in Table [H represent a compromise between trap sta- 
bility and anharmonic effects The Fermi temper- 
ature is given by Tp = Ep/k, where the Fermi energy 
Ep = h{^NujxUJyUJzy^^ = h'^kp/2m, kp is the Fermi 
wavenumber and k is the Boltzmann constant. The pa- 



Table I: Trap parameters for the different modes. 
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Figure 2: Sloshing mode frequency uos normalized by the trap 
frequency uJx as a function of temperature. The measured fre- 
quency shows a decrease with increasing temperature, (dots) 
which is due to the increase in the size of the cloud. The lines 
show the expected frequency from a first-order anharmonic 
correction; see Appendix B. To determine the cloud size for 
different temperatures we assume a harmonic potential (solid 
line) and a Gaussian potential (dashed line) , respectively [l^] ■ 



rameter Vq is the trap depth, and N is the total number of 
atoms, given by iV = 6 x 10^. The interactions are char- 
acterized by the dimensionless parameter 1/kpa, where 
a is the s-wave scattering length. 

To control the aspect ratio uj^jiOy^ we use rapid spatial 
modulation of the trapping beam by two acousto-optical 
deflectors, resulting in the creation of time-averaged trap- 
ping potentials This on one hand allows us to com- 
pensate for residual ellipticity of the trapping potential 
on the percent level and thus to realize a cylindrical sym- 
metric trap = Wy). On the other hand it allows 
for the excitation of surface modes by deliberately in- 
troducing elliptic trapping potentials [i^x 7^ ^y)- The 
procedures used to excite the modes are outlined in Ap- 
pendix A. To change the temperature we apply a con- 
trolled heating scheme via sudden compression of the gas 
as described in 10]. Detection of the cloud is done by 
absorption imaging which displays the shape of the cloud 
in the x-y plane after expansion. For each mode under 
investigation we determine the frequency and damping 
following the procedures of our previous work lol. llOj] ; 
see also Appendix A. 

Because of the Gaussian shape of the trapping poten- 
tial, corrections are needed for a precise comparison of 
the experimental observation to the idealized case of per- 
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feet harmonic trapping. Especially for higher tempera- 
tures, when the size of the cloud is larger, anharmonic 
corrections become important. This is demonstrated 
by measurements of the transverse sloshing mode fre- 
quency ujs (Fig. [21) , which clearly show a substantial de- 
crease with increasing temperature. To reduce the anhar- 
monic effects on the frequencies of the collective modes 
under investigation, we normalize the compression and 
quadrupole collective mode frequencies to the sloshing 
mode frequency in the transverse direction. This nor- 
malization reduces the anharmonic effects to a large ex- 
tent since the decrease of the sloshing mode frequency 
with increasing cloud size is of the same order as the cor- 
responding decrease of the frequency of the transverse 
modes jl6l | . To normalize the scissors mode frequency we 
take the geometric average of the two different sloshing 
mode frequencies in the transverse direction. 

For each of the trap parameters of the different modes 
we determine the sloshing mode frequency as a function 
of the temperature. As an example we show lOs for the 
trap parameters used for the compression mode measure- 
ment; see Fig. [21 We compare ujs/^x (dots) to a theoreti- 
cal model which allows to calculate the sloshing frequency 
as a function of the cloud size; see Appendix B. Assuming 
a harmonic potential to derive the mean squared size (a;^) 
[iBl underestimates the anharmonic effects (solid line) in 
particular for higher temperatures. Taking into account a 
Gaussian potential to determine {x'^) (dashed line) agrees 
much better with the measured sloshing frequency. 



III. THEORY 

We shall compare our experimental findings to the re- 
sults of model calculations that apply to the normal state 
of the gas, i.e. at temperatures above Tc- In this Sec- 
tion, we outline our theoretical approach to the calcula- 
tion of mode frequencies for T > Tc- A more detailed 
description can be found in Refs. [13] and 21 1. We as- 
sume that single-particle excitations are reasonably well 
defined in the sense that most of the spectral weight of 
the single-particle spectral function is located at a peak 
corresponding to that of non-interacting particles. The 
low-energy dynamics of the gas can then be described 
by a semiclassical distribution function f{r,p,t) which 
satisfies the Boltzmann equation. A collective mode cor- 
responds to a deviation Sf = f — f'^ away from the 
equilibrium distribution /'•\r,p). Writing Sf{r,p,t) = 
/"(r, p)[l — /'^(r, p)]<i>(r, p, t) and linearizing the Boltz- 
mann equation in 6f(r,p,t) yields 



rO^ fd^ -9^ -9^ 



-mh (1) 



where r = v = p/m, p — —dV/dr and / is the collision 
integral. We take the potential V{r) to be harmonic and 
given by V{r) = ■m{uj1x'^ + J^y^ -I- uJ\z^)l2. 

To describe the collective modes we expand the devi- 
ation function in a set of basis functions according 



to 



$(r,p,i) = e 



«(r,p), 



(2) 



Since the purpose of this article is the comparative 
study of different collective modes and not the precision 
measurement of a single mode as in previous work , we 
follow a faster yet simpler procedure to normalize the fre- 
quencies. We measure the sloshing mode frequency only 
at particular temperatures of interest. From these points 
we determine the temperature dependence of the slosh- 
ing frequency by interpolation. Even though the normal- 
ization takes into account the temperature dependence of 
the anharmonicity, it does not reduce effects due to drifts 
in the power of the trapping beam. We believe this to be 
the main source for the scatter of the data in Fig. [4l 

To determine the temperature of the gas we first adia- 
batically change the magnetic field to 1132 G [l3], where 
l/kpa ~ —1, to reduce the effect of interactions on 
the density distribution [l^. Under this condition, for 
T > 0.2Tf, the interaction effect on the density dis- 
tribution is sufficiently weak to treat the gas as a non- 
interacting one to determine the temperature from time- 
of-flight images. We fit the density distribution after 2 ms 
release from the trap to a finite-temperature Thomas- 
Fermi profile. The temperature measured at 1132 G is 
converted to the temperature in the unitarity regime un- 
der the assumption that the conversion takes place isen- 
tropically, following the approach of Ref . [T^ . Statistical 
uncertainties for the temperature stay well below O.ObTp. 



where w is the mode frequency. For the compression 
mode with a velocity field v cx {x,y,cz), with c a con- 
stant, we use the functions 

= + y^, (1)2 = xpx + ypy, (1)3 = pI+pI, (1)4 = pI- 

(3) 

For the quadrupole mode with a velocity field v cx 
(x, — y, 0) (ignoring the small velocity along the axial di- 
rection), we use 

01 = - y^, (1)2 = xp^ - ypy, (f)s=pl- pI, (4) 

whereas the basis functions for the scissors mode are 
given in Ref. [l^l . Our choice of basis functions is physi- 
cally motivated by the characteristic features of the three 
different modes illustrated in Fig. [1] Since we limit our- 
selves to a few simple functions, the basis sets are not 
complete, but we do not expect qualitative changes to 
occur as a result of including more basis functions in our 
calculation. 

We now insert the expansion ^ into ^ and take 
moments by multiplying with the functions (j)i and in- 
tegrating over both r and p. This yields a set of linear 
equations for the coefficients q for each of the collective 
modes. The corresponding determinants give the mode 
frequencies. For the compression mode, we obtain 
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Figure 3: The effective collision rate for a gas in the unitarity 
limit. The dashed curve is the classical result, the dash-dotted 
includes Pauli blocking, and the solid line includes pairing 
correlations in the scattering matrix. The superfluid region 
for T < Tc is indicated. 



and for the quadrupole mode, we get 



(6) 



The equation for the scissors mode is given in Ref. |23 |. 
The effective colhsion rate l/r in ([5]) and jl]) is given 

by 



/ d^rd^ppxPyI[pxPy] 

jd?rcpppiplp{\^ n 



(7) 



Note that this expression for 1/t involves a spatial av- 
erage over the cloud. In the coUisionlcss Umit, ujt ^ 1, 
the two equations ([5]) and © both yield u; = 2uj±, where 
= oJx = while in the hydrodynamic limit, ldt <^ 1, 
they result in w = v/l0/3wi_ for the compression mode 
and Lu — ^/2uj± for the quadrupole mode. 

The dependence on temperature T and scattering 
length a enters through t. In particular, Pauli blocking 
and pair correlations strongly depend on T and a, and we 
now examine their role on the effective collision rate. In 
Fig. [3l we plot 1/t as a function of temperature for a gas 
in the unitarity limit |a| ^ cxd using three different ap- 
proximations for the collision integral. First, the dashed 
curve gives the effective collision rate in the classical 
regime using the vacuum expression %a,c = ^/(l + iqa) 
for the scattering matrix with Tq = 47rS^a/m. The s- 
wave differential cross section da/d^l which enters in the 
collision integral / is related to the scattering T matrix 
by da/dQ = m'^\T\'^ / (ATrh'^f . In the classical regime, we 
then get from ([7]) 
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4 kTpTj, 

Tclass 457r h T2 



(8) 



for a gas in the unitarity limit j22j . The small prefactor 
4/(457r) « 0.028 in ^ imphes that the effective colhsion 



rate is significantly smaller than what one would expect 
from simple estimates or dimensional analysis at unitar- 
ity. Second, the dash-dotted curve gives the effective 
collision rate when Pauli blocking effects are included as 
in [20] , while we still use the vacuum expression T^ac for 
the scattering matrix. Pauli blocking effects reduce the 
available phase space for scattering thereby reducing the 
scattering rate. For T <C Tp Pauli blocking in a normal 
Fermi system gives l/r oc T^. Finally, we plot as a solid 
curve in Fig. [3] the effective collision rate taking into ac- 
count both Pauli blocking and many-body effects for T 
in the ladder approximation which includes the Cooper 
(pairing) instability. This gives T — %/{! — TqU) where 
n is the pair propagator. Since our treatment of the pair 
correlations only apply to the normal state of the gas, we 
plot this curve for temperatures greater than the critical 
temperature Tc, which within the ladder approximation 
used here is given by Tc « O.STp for a trap f2ll|. 

We see that 1/t is increased by the pairing correla- 
tions for the T-matrix. The pairing correlations signif- 
icantly increase the effective collision rate for tempera- 
tures (T - Tc)/Tc < 1 HH. One often refers to this tem- 
perature range as the pseudogap regime. In fact, pair- 
ing correlations almost cancel the Pauli blocking effect in 
the collision integral above Tc and 1 /r is fairly accurately 
given by the classical value as can be seen from Fig. [3l At 
high temperatures, this cancelation can be demonstrated 
analytically by carrying out a high-temperature expan- 
sion of ([7]). We obtain after some algebra the simple 
expression 



1 



'^class 



l_ (Tf 
32 V T 



(9) 



The presence of the small prefactor 1/32 in ^ shows 
that the leading correction to the classical limit is less 
than 3% at temperatures above the Fermi temperature 
Tp. 



IV. RESULTS AND DISCUSSION 

The theoretical results of the previous section were 
all obtained for a purely harmonic potential. Since an- 
harmonicity plays an important role in our experiments, 
as discussed in Sec. |TT1 we normalize the measured fre- 
quencies and damping rates for the collective modes to 
the measured temperature-dependent sloshing frequen- 
cies, for which an example is shown in Fig. O In the 
following we compare our observations to the theoretical 
results. It should be emphasized that the theoretical ex- 
pressions for the frequency and damping contain no free 
parameters to fit theory and experiment. 

First we discuss the frequency for the three modes un- 
der investigation as a function of the temperature, as 
plotted in Fig. 21 In all three cases the theoretical expres- 
sion for the frequency (the full lines in Fig. 2]) smoothly 
changes from the hydrodynamic value at the lowest tem- 
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Figure 4: The three panels show the observed normal- 
ized mode frequencies versus temperature for the quadrupole 
mode, the scissors mode and the compression mode. The 
error bars indicate the statistical error of a single frequency 
measurement. The full lines are the result of the theory for 
the normal state described in the section Hill which includes 
the combined effects of Pauli blocking and pair correlations; 
note that these curves start at T = O.STf, which in the ladder 
approximation used here is the transition temperature to the 
superfluid state. For illustrative purposes we also show the 
theoretical results when only Pauli blocking is taken into ac- 
count (dash-dotted lines) and those for a classical gas (dashed 
lines) . 



perature considered to the coUisionlcss value at high tem- 
peratures. The normahzed frequencies in the hydrody- 
namic limit for the quadrupole mode and compression 
mode are -\/2 « 1.41 and y'lO/S 1.83, respectively. 
The normalized frequency in the coUisionless hmit for 
both these modes is 2. Using the geometric average of 
the trap frequencies to normahze the scissors mode fre- 
quency, one gets, using the ratio ujx/ujy = 16/7 from 

Table I, that ' 



- ujy) / (ujx^y) ~ 1-65 in the hydro- 
dynamic hmit and {ujx + ujy) / ^uixOJy ~ 2.17 in the coUi- 
sionless limit. Note that the scissors mode consists of a 
two- frequency oscillation in the coUisionless limit. Here 
we only consider the larger frequency component. The 
lower frequency component exhibits increasing damping 
towards lower temperatures and disappears in the hydro- 
dynamic limit j23 |. 

Figure [5] illustrates that there is a reasonable over- 
all agreement between experiment and theory, although 
some differences exist. The agreement is best for the scis- 
sors mode, while for the quadrupole mode the changeover 
from hydrodynamic to coUisionless behavior happens at a 
lower temperature than the one found theoretically. The 



03 

"S 

O) 

c 

Q. 

E 

"O 

T3 
(U 
N 

TO 

E 
\— 
o 
c 



0.4 

0.2 

0.0 
0.4 

0.2 

0.0 
0.10 

0.05 



0.00- 




quadrupole mode 4 

• y • • 




scissors mode 




compression mode 





0.2 0.4 0.6 0.8 1.0 1.2 

T/T. 



Figure 5: Normalized mode damping versus temperature for 
the quadrupole mode, the scissors mode and the compres- 
sion mode. The points are experimental values, while the full 
lines represent our calculated values, taking into account both 
Pauli blocking and pairing effects. The dash-dotted line only 
takes Pauli blocking into account, while the dashed line is the 
classical (high-temperature) result. 



measured compression mode frequency, which shows con- 
siderable scatter, increases with increasing temperature 
and is close to the coUisionless value at the highest tem- 
perature measured. 

The observed change from the hydrodynamic to the 
coUisionless frequency for the compression mode is in 
contrast to Ref . [7| , where the frequency remains close to 
the hydrodynamic value for the same temperature range. 
We attribute this discrepancy to different treatments of 
anharmonic effects, which are particularly important for 
this mode since the difference between the hydrodynamic 
and coUisionless frequency is of the same order as the 
frequency shift due to anharmonic effects. In Ref. [7*1 the 
data are corrected by including anharmonic effects to first 
order, while we adopt the point of view that the main an- 
harmonic effects can be taken into account by normaliz- 
ing the measured oscillation frequencies to the measured 
temperature-dependent sloshing frequencies. Fig.[2]iUus- 
trates that a simple first-order treatment of anharmonic 
effects on the sloshing frequency does not account quan- 
titatively for the observed variation with temperature. 

At very low temperatures the measured frequencies are 
close to the hydrodynamic values because the gas is in 
the superfluid phase [1, [^. Without pair correlations, 
but with Pauli blocking, at these low temperatures the 
frequencies would assume their coUisionless values as il- 
lustrated by the dashed-dotted lines in Fig. [H 

We now proceed to consider the damping of the os- 
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cillations. The experimental values for the normalized 
damping rate are shown in Fig. [5l Theoretically, one 
expects the damping to vanish in the hydrodynamic and 
collisionless limits and exhibit a maximum in between, as 
brought out by the calculations in Sec. IIIII Experimen- 
tally, both the quadrupole and the scissors mode exhibit 
the expected maximum in damping in the transition re- 
gion. For the compression mode, however, the damping 
does not decrease at higher temperatures. This surpris- 
ing behavior for the compression mode has already been 
found in 0. A possible reason for the increasing damp- 
ing is dephasing-induced damping due to anharmonicity. 
Anharmonic effects are more important for the compres- 
sion mode as the intrinsic damping is relatively small 
due to the small difference between the frequencies in 
the collisionless and hydrodynamic limits [2J|. In con- 
trast to the case of frequency discussed above, we can- 
not expect to take into account the main effects of an- 
harmonicity by normalizing the measured damping rate 
to the temperature-dependent sloshing mode frequencies. 
This makes it delicate to compare our experimental re- 
sults to those of a theory based on a purely harmonic 
potential. The damping of the quadrupole mode shows 
the expected qualitative behavior, although the maxi- 
mum in damping happens at a lower temperature com- 
pared to theory. This is consistent with the frequency 
data for this mode, since the transition there also hap- 
pens at lower temperature. For the scissors mode the 
experimental data agree fairly well with theory, although 
some discrepancy exists at the lowest temperatures. 

We can relate the frequency and damping of the 
quadrupole mode directly to each other by eliminating 
the collision rate l/r in ([5]) . Writing lo = ojq — iTq for 
the solution of ^ with luq and Tq being the quadrupole 
frequency and damping, we obtain 
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Figure 6: Normalized mode damping versus normalized fre- 
quency for the quadrupole and scissors mode. The solid line 
shows the expected behavior for a harmonic trap. The arrows 
point toward the direction of increasing temperature. 



V. CONCLUSION 



Fq = J-^l - ^1 + ^Su^l - 7u;l. (10) 



A similar relation holds for the two other modes. This 
allows us to compare theory and experiment indepen- 
dently of any approximations involved in the evaluation 
of l/r. Figure [6] shows the normalized damping rate ver- 
sus the normalized frequency of the quadrupole and the 
scissors mode; we do not show the data for the com- 
pression mode because of the apparent problems dis- 
cussed before. We find that the maximum damping of 
the quadrupole mode is larger than expected. For the 
scissors mode the damping is larger only at low frequen- 
cies. This suggests that the difference between theory 
and experiment is not a consequence of the approxima- 
tions entering the calculation of the relaxation rate but 
could be due to anharmonic effects or the need for larger 
basis sets to describe the modes [see and (|H)]. 



In this work we have presented measurements of the 
frequency and damping of three different collective modes 
under similar conditions for an ultracold Fermi gas of ^Li 
atoms in the unitarity limit. The experimental results 
obtained in the normal state of the gas are in reason- 
able agreement with our theoretical calculations, which 
take into account Pauli blocking and pair correlations. 
The remaining discrepancies may originate in a variety 
of sources such as our treatment of anharmonic effects, 
the temperature calibration, and the use of a restricted 
basis for solving the Boltzmann equation. Also they may 
reflect the need to incorporate further interaction effects 
in the kinetic equation, which forms the starting point 
for the theoretical calculations. For instance, there are 
self-energy shifts on the left-hand side of the kinetic equa- 
tion which could be important. The study of collective 
modes is a sensitive probe of the properties of strongly in- 
teracting particles such as the gas of ^Li atoms under in- 
vestigation, and further work on temperature-dependent 
phenomena will undoubtedly shed more light on these 
interesting many-body systems. 
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Appendix A 

Here we present more details on the experimental pro- 
cedures to excite the three collective modes. 

To excite the radial quadrupole mode we adiabatically 
deform the radially symmetric trap to an elliptic shape 
while keeping the average trap frequency constant be- 
fore turning off the deformation suddenly The de- 
formation is chosen such that the amplitude of the mode 
oscillation relative to the cloud size is below 10%. A two- 
dimensional Thomas-Fermi profile is fitted to the images, 
taken after a short expansion time of 0.5 ms. The dif- 
ference in the width of the main axes is determined for 
different hold times and fitted to a damped sine function, 
from which we determine the frequency and damping of 
the mode. 

The excitation of the radial compression mode is done 
by a sudden compression of the cloud. To determine 
the frequency and damping of the compression mode we 
follow the same procedure as for the quadrupole mode 
but fitting to the sum of the widths. Here we use an 
expansion time of 2 ms before taking the image. 

The scissors mode appears as an angular oscillation of 
an elliptic cloud about a principal axis of an elliptic trap. 
To excite this oscillation we create an elliptic trap in the 
x-y plain and suddenly rotate the angle of the principal 
axes by 5 degrees [l3|- The tilt of the principal axes of 



the cloud is determined 0.8 ms after releasing the cloud 
from the trap for different hold times. If the gas is hy- 
drodynamic, we fit a single damped sine function to the 
oscillation of the angle. However, for a coUisionless gas, 
the oscillation exhibits two frequencies. Thus wc fit a 
sum of two damped sine functions each with their own 
free parameters. When the behavior changes from hydro- 
dynamic to coUisionless the single damped sine function 
fits the data reasonably well, as discussed in Since 
the larger of the two frequencies in the coUisionless regime 
smoothly connects to the hydrodynamic frequency at low 
temperatures we only consider this frequency in the pa- 
per. 



Appendix B 

Here we briefly discuss the calculation of the transverse 
sloshing modes including anharmonic corrections to low- 
est order. The transverse trapping potential is 

(x^ x'^y^ \ 

^+6^"2^"2&^"^j- (^^^ 

Concentrating without loss of generality on the slosh- 
ing mode in the x-direction, we choose the function 
— cix + C2Px- Putting this into the linearized Boltz- 
mann equation ([1]), eliminating C2, and taking the mo- 
ment J dxdyn{x, y) with n[x,y) the density (we ignore 
the axial direction) , we obtain for the sloshing frequency 

Here {x^) = J n{x,y)x'^dxdy / J n{x,y)dxdy and we have 
used Lol — 2VQ/ma? together with ujy = 2VQ/mb'^. 
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